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Demagnetization, which is inherently present in the magnetic response of small finite-size super- 
conductors, can be accounted for by an effective k within a two-dimensional lowest Landau level 
approximation of the Ginzburg-Landau functional. We show this by comparing the equilibrium 
magnetization of superconducting mesoscopic disks obtained from the numerical solution of the 
three-dimensional Ginzburg-Landau equations with that obtained in the "effective" LLL approxi- 
mation. 



I. INTRODUCTION 

The magnetic response of a three-dimensional super- 
conductor depends on intrinsic parameters such as the 
coherence length ^ and the penetration length A through 
the ratio k — A/^, but also depends critically on the 
shape and orientation with respect to the applied mag- 
netic field H . The latter dependence makes it very diffi- 
cult to extract information about the true nature of the 
superconductor from magnetic measurements and it can 
completely mask the intrinsic magnetic response. This 
phenomenon is known as demagnetization. The best 
known example of how demagnetization affects the mag- 
netic properties of a finite-size superconductor is the ap- 
pearance of the intermediate state |l| in type-I super- 
conductors. In very simple terms the intermediate state 
developes due to the tendency for the expelled and "over- 
stretched" magnetic flux lines to penetrate back into the 
otherwise perfectly diamagnetic material. For instance, 
macroscopic type-I superconducting disks in perpendicu- 
lar geometry or finite-length cylinders with the rotational 
axis parallel to the field, which are still the focus of exper- 
imental and theoretical studies [^[j4j , exhibit an apparent 
type-II magnetic response due to the formation of the in- 
termediate state. Thin mesoscopic Aluminum disks (with 
radii ranging typically between 0.1 and 2/Ltm) have also 
received a lot of attention lately, in part due to recent 
breakthroughs in magnetic |^ and resistive |^ measure- 
ment techniques. Demagnetization in these mesoscopic 
disks is the focus of this paper. 

From the theoretical point of view the demagnetization 
presented by macroscopic ellipsoidal samples can be sim- 
ply accounted for by a demagnetizing factor, N , which 
is a scalar if H lies parallel to one of the principal axis 
of the ellipsoid [Q. N runs from N = 1 (maximum de- 
magnetization) for the case of an infinite ellipsoid per- 
pendicular to to TV = (zero demagnetization) for 



the same ellipsoid parallel to H , taking any value in be- 
tween for intermediate situations like the spherical ge- 
ometry {N = 1/3). When the samples are not ellipsoidal 
the situation is more complicated. Still, various effective 
models have been proposed to account for demagnetiz- 
ing effects in macroscopic samples with simple forms like 
strips, disks or finite cylinders with various orientations 
with respect to the field [^] , and very accurate analyti- 
cal results have been obtained in certain limits jj] . Com- 
pared to macroscopic disks, the situation for mesoscopic 
disks 1^ gets complicated by the fact that A can be com- 
parable to the dimensions of the sample; in particular, to 
the dimension parallel to the field or thickness of the disk, 
d. When this is the case the penetration length is naively 
expected to get renormalized to A = A^/d and the new 
K can be larger than l/\/2, value that separates type-I 
from type-II behavior. Whether or not thin mesoscopic 
Aluminum disks exhibit truely type-II behavior and vor- 
tices form in the condensate can oiJy be unambiguously 
answered with imaging techniques ||8|| since demagnetiza- 
tion is always present and complicates the interpretation 
of the experimental results ||]. 

In principle, this non-trivial interplay between di- 
mensions, geometry, and intrisic parameters can only 
be addressed theoretically within a three-dimensional 
Ginzburg-Landau framework: 

_1_ (^.,nV - ^^(r)) *(f) = -a^{r) - [im{r)\m{r)\\ 

(1) 

^ ^ ^ 47r^ 
V X V X A{r) = — j, (2) 
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where j is the superconducting current |l| . These equa- 
tions are obtained after extremizing the phenomenologi- 
cal Ginzburg-Landau energy functional: 
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between the exact results and those obtained in our "ef- 
fective" LLL approach. 
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where Gn is the Gibbs free energy of the normal state and 
—ihSI — e* A(f)lc is the momentum operator for Cooper 
pairs of charge e* = 2e and mass m* = 2m in a vec- 
tor potential A(f) which is associated with a magnetic 
induction h,(f). is the Cooper pair wave function 

or complex order parameter and the coefficients a and /? 
have the usual values [Q which scale phenomcnologically 
with temperature . 

The difficult task of solving the coupled Ginzburg- 
Landau differential equations (|l|) and (||) for three- 
dimensional mesoscopic disks has been undertaken nu- 
merically by Peeters and co-workers (see also Ref. | |To[ ] 
for the equivalent problem of cylinders). While, typically, 
the order parameter within the disk can be considered 
uniform in the direction of the field, i.e., effectively two- 
dimensional, the three-dimensional magnetic flux struc- 
ture needs to be taken fully into account. Consider- 
ing the expected limitations of the effective Ginzburg- 
Landau theory, the resultsjp| reproduce to a large extent 
the experimental findings pT Akkermans and co-workers 
have undertaken similar studies from a different point 
of view. They have been able to minimize analytically 
the Ginzburg-Landau functional close to the dua.1 point 
K ~ and in the London limit k — > oo [|l2| , al- 

though demagnetizing effects have been considered only 
at a phenomenological level. 

The aim of this work is to show that for small disks 
there is an effective way to incorporate demagnetization 
in the solution of the Ginzburg-Landau equations with- 
out considering in detail the three dimensional structure 
of the magnetic field lines. Essentially, we do this by 
introducing an effective k, which we will denote by k, 
that is geometry dependent. Adequately chosen, the use 
of K instead of the intrinsic n can account, in part, for 
demagnetizing effects. This is conveniently done within 
the framework of the lowest Landau level (LLL) approx- 
imation which has been broadly used in the case of bulk 
superconductors ]l^ , p^ and mesoscopic systems p5| , p^ . 
It will turn out that this procedure works rather well if 
we define two different effective k, one for the low mag- 
netic field region, i.e., the Meissner state, and one for 
magnetic fields near Hc2- 

In Sec. H we briefly review how the standard LLL ap- 
proximation must be modified to obtain the magnetic re- 
ponse of finite-size superconductors where surface effects 
are dominant. We also compare our LLL approximation 
with traditional JlTf results for magnetization in the case 
of bulk systems. In Sec. Ill we present results for the 



equilibrium magnetization of mesoscopic disks obtained 
from the full numerical solution of the Ginzburg-Landau 
differential equations [^. We finally introduce the effec- 
tive K in our LLL approximation and make a comparison 



II. MAGNETIZATION IN THE LOWEST 
LANDAU LEVEL APPROXIMATION 

A. The basics 

A fully two-dimensional alternative to the approaches 
mentioned in the introduction for solving the Ginzburg- 
Landau equations has been proposed by one of the au- 
thors [ p^ . One expands the order parameter in a set of 
functions that are the lowest level solutions of the lin- 
earized version of Eq. (pi) : 
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It is essential for the radial part of the expansion func- 
tions, $L(r), to obey the standard boundary conditions 
of zero current through the surface The presence 

of the boundary implies that, for large enough angular 
momentum L, an infinite number of bulk Landau levels 
participate (numerical work is usually required at this 
point). Nevertheless, we will still refer to this expansion 
as a LLL expansion. The other central idea behind this 
approximation is to consider the magnetic induction to 
be spatially constant, B, but not necessarily equal to the 
external field H , as it has been usually the case in the 
context of the LLL approximation for bulk systems [ p^ . 
The magnetic inducti on B sets the scale of the magnetic 
length £ — y/e*hJcB and, thus, the scale of the expan- 
sion functions. Obviously, to consider a constant value 
of B across the superconductor is strictly valid only in 
the limit B H and it does not capture correctly the 
magnetic induction profile in the Meissner phase. On the 
other hand, the LLL approximation is not expected to be 
valid below B « Q.3Hc2, although this naive expectation 
has been shown to be higher than the real limit . Still, 
as we will show in the next section, this approach gives 
qualitative and quantitatively good results for small and 
thin Aluminum disks in the whole range of fields as long 
as K is adequately chosen. 

When one substitutes this expansion into the 
Ginzburg-Landau functional (|^) one obtains an algebraic 
expression in terms of the complex coefficients Cl- 
og 

Gs = Gn + Y.{a + eL)\CL\'' 

-0 

OO „ 

J2 CICI^Gl.Cl, / dfTlJljL.TL, 

_ r _ T n 



L=0 



2 



Z/i,L2,-/v3,L4 — 

'V{B ^ Hf /St:. 



(5) 



The effect of the surface on the condensate is cast in 
the Cooper pair energy e^, which, unlike bulk systems. 
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presents a dip for states L close to the surface [|l5[, and 
in the "interaction" integrals. In this form the Ginzburg- 
Landau functional can be easily minimized or, in general, 
extremized with respect to C^. and B. This can even 
be done analytically when not too many terms in the 
expansion (Q) participate 15 1^. 



B. Bulk magnetization 

It is illustrative at this point to consider a bulk sys- 
tem within our LLL approximation. Let us assume an 
ellipsoidal form or a long cylinder for which no demag- 
netization exists {N = 0). In the spirit of Abrikosov's 
seminal work, it is straightforward to obtain an expres- 
sion for the magnetic induction 
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and the free energy density 
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In the above expressions V is the volume and we 
have made use of the Abrikosov parameter (3a = 
(^'(r)2)'^/(^'(r)^) |l^ which measures the uniformity of 
the superconducting density. Notice that the expression 
that we obtain for the magnetic induction is different 
from that obtained by Abrikosov [pT[: 
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The difference lies essentially in two facts. (i) We 
consider a uniform magnetic induction B while, in 
Abrikosov's work, h(f) is not uniform; instead, the spa- 
tial average (/i(r)) is calculated, (ii) Our LLL expan- 
sion functions are evaluated for a value of B that needs 
to be determined after minimization; on the other hand, 
Abrikosov's expansion functions are calculated at the up- 
per critical field Hc2- These two significant differences ac- 
count for the difference in the magnetization expressions 
although, our approach is considerably simpler. Consid- 
ering that when one minimizes with respect to the struc- 
ture of the vortex lattice one obtains Pa ~ 1.16, which 
corresponds to a triangular arrangement of vortices [Q, 
both approaches give a remarkably similar result for the 
bulk magnetization. 



III. DEMAGNETIZATION IN THE LOWEST 
LANDAU LEVEL APPROXIMATION 

We now compare our LLL approximation with the ex- 
act magnetization results for mesoscopic disks of typical 
sizes i? ~ ^. Here, only one state T/, participates in the 
expansion (0) for any value of the external field or, in 



other words, the order parameter has a well-defined value 
of the angular momentum and forms a giant vortex [po[ . 
Minimizing Eq. ^ with respect to Cl one obtains 
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where the minimal value of B is obtained numerically. 
The free energy is given in units of iJ^jl^/STr, eL{B) 
is given in units of the LLL energy huJc/2 (with lOc = 
e*B/m*c), R is in units of ^, and B and H are given in 
units of Hc2- Notice that in the denominator of the above 
expression we have written k instead of k. As mentioned 
in the introduction, finite thickness affects the real value 
of K through the expression A^/d, but also does demag- 
netization in a more complicated way. We propose to 
consider k as an effective parameter that depends on the 
geometry, i.e., on R and d in such a way that the exact 
magnetic response is approximately reproduced (within 
the limits of the LLL approximation). Figure |^ shows 
the equilibrium pT|] magnetization of a disk character- 
ized by i? = 3^ and k = 1 for four different values of 
the disk thickness: d/£_ = 0.01(a), 0.1(6), 0.5(c), 1.0(d). 
Dashed lines correspond to the LLL approximation and 
solid lines to the three-dimensional numerical solution of 
Eqs. @) and d). 
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FIG. 1. Equilibrium magnetization of a disk of radius 
i? = 3^ and bulk k = 1 for four different thicknesses; 
dji = 0.01(a), 0.1(6), 0.5(c), 1.0(d). Solid lines correspond to 
the full numerical solution and dashed lines correspond to the 
LLL approximation using different values of k: (a) k = 32, 
(b) k = 6.1, (c) k = 2.9, and (d) k = 2.15. 

Each jump is associated with a transition from one 
giant vortex to another with one unit more or less of 
vorticity (L ^ L ± 1) ||l|,|,@. The value of k has been 
chosen to match the slope of the quasi-linear magnetic 
response of the Meissner state (L = 0). As mentioned 
before, it is in this state that our approximations are 
expected to be less accurate since the magnetic induction 
is clearly not uniform there; still, keeping in mind that 
we have chosen the worst case scenario for our fitting. 
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it is a simple and unambiguous criterion for choosing k 
and we will use it in what follows. In this figure one 
can appreciate that the magnetization obtained in the 
LLL approximation follows closely the "exact" magnetic 
response, both in magnitude and position of the jumps. 
As the thickness increases and the disk resembles more 
and more a long cylinder, the value of k approaches the 
intrinsic one and the effective LLL approximation works 
better. Notice, however, how demagnetization manifests 
itself in that the value of k [(a) k = 32, (b) k = 6.1, 
(c) k = 2.9, and (d) k — 2.15] is different from the one 
expected from the simple scaling A^/d (see Fig. ||). This 
manifestation is more clear for larger disks as will be 
shown below. 




FIG. 2. Equilibrium magnetization of a disk of thick- 
ness d = 0.1^ and bulk k = 0.28 for different radii: 
= 2.0(a), 3.0(6), 4.0(c), 5.0(d). Solid lines correspond to 
the full numerical solution, dashed lines correspond to the 
LLL approximation using different values of k: (a) k = 1.45, 
(b) k = 1.78, (c) k = 2.2, and (d) k = 2.7, and dotted lines 
correspond to the LLL approximation using (a) k = 1.20, (b) 
k — 1.35, (c) k = 1.45, and (d) k = 1.55 

Confident in the validity of the effective LLL approxi- 
mation to describe the actual magnetic response of small 
type-II disks (we have repeated the analysis just de- 
scribed for different disk radii with similar encourag- 
ing results), we now try to push it further to include 
the description of type-I mesoscopic disks like those of 
Geim's experiments . Figure shows the equilibrium 
magnetic response of an Aluminum disk (k — 0.28) of 
thickness d = 0.1^ for four different values of R/^ ~ 
2.0(a), 3.0(6), 4.0(c), 5. 0((i). The same convention as be- 
fore for fixing k has been followed (dashed lines). Al- 
though for d = 0.1^ a renormalized k is already expected, 
it is clear from the results that, as the radius increases, 
one needs to increase k [(a) k = 1.45, (b) k = 1.78, (c) 
k = 2.2, and (d) k — 2.7] in order to match the posi- 
tion of the magnetization jumps and, to some extent, the 
magnitude of them. The fact that we have to use higher 
values of k as the disk radius increases (whereas the thick- 



ness remains constant) constitutes again a manifestation 
of demagnetizating effects which decrease the magnetic 
response of the disk as the radius grows bigger (i.e., the 
disk resembles more and more a thin film in perpendic- 
ular geometry). Even for the larger disk considered the 
agreement of the LLL results with the exact solution is 
fairly satisfactory. However, now it becomes more diffi- 
cult to match the overall magnitude of the magnetization 
in the whole range of fields with a single value of k. At 
large fields, where the LLL approximation is expected to 
work better, the agreement can be improved considering 
a different set of values for k (dotted lines). Unlike the 
previous case, the fitting criterion is a little ambiguous at 
high fields, but one can clearly appreciate the improve- 
ment. As Fig. ^ shows, in all cases considered k behaves 
almost linearly with R. One last caveat: Multiple- vortex 
states could be stable in the condensate of the larger 
disk, but we are not considering here such a possibility 
since it is not relevant for our discussion. 
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FIG. 3. Equilibrium magnetization of a disk of thickness 
d = 0.5^ and bulk k = 0.28 for different radii. Solid lines 
correspond to the full numerical solution, dashed lines corre- 
spond to the LLL approximation using different values of k: 
(a) k = 0.75, (b) k = 0.95, (c) k = 1.18, and (d) k = 1.45, 
and dotted lines correspond to the LLL approximation using 
(a) k = 0.55, (b) k = 0.65, (c) k = 0.70, and (d) k = 0.75 

Finally, Fig. ^ shows the equilibrium magnetization 
of an Aluminum disk with d — 0.5^ for different radii: 
R( = 2.0(a), 3.0(6), 4.0(c), 5.0(d). Now it becomes more 
difficult to find an effective k that allows the LLL ap- 
proximation to reproduce, even if only qualitatively, the 
exact results. This is not surprising since, for such thick- 
ness, the type-I behavior of Aluminum manifests itself 
more pronouncely. Still, the number of magnetization 
jumps is approximately reproduced in all cases shown, 
but the the magnitude of the magnetization is apprecia- 
bly underestimated with the set of values for k that fit the 
Meissner phase (dashed lines). Again, a different set can 
be chosen so that the high-field response is approximately 
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reproduced (dotted lines), although the agreement is not 
completely satisfactory for large disks. 

As conclusion, and for illustration purposes, we plot 
in Fig. 4 function of thickness and radius, in- 

cluding the cases considered above. With the help of the 
calculations presented in this paper, we leave it to the 
reader to judge by himself when the effective LLL can be 
safely used and when and how demagnetization effects 
can be accounted for in this simple manner. It is beyond 
the scope of this paper to present a thourough study of 
all possible cases, but we expect this work to serve as a 
useful guide to the specialized reader. 




FIG. 4. (a) Black dots: function of the thickness for 

n — 1 and R — 3^. White dots: k expected from the scaling 
X^/d. (b) « as a function of R for k = 0.28 and d — 0.1^ 
(upper curves) and d = 0.5^ (lower curves). Black symbols 
denote fitting to the Meissner phase and white symbols to the 
high-field magnetic response. 
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